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Abstract 

In this paper we address practical aspects of the implementation of the 
0-1 test for chaos in deterministic systems. In addition, we present a new 
formulation of the test which significantly increases its sensitivity. The test can 
be viewed as a method to distill a binary quantity from the power spectrum. 
The implementation is guided by recent results from the theoretical justification 
of the test as well as by exploring better statistical methods to determine the 
binary quantities. We give several examples to illustrate the improvement. 



1 Introduction 

Being able to distinguish between regular and chaotic dynamics in a deterministic 
system is an important question with applications ranging from cardiac arrhythmias 
to the stability of our solar system. Much progress has been made in developing 
tests for chaos [13], HH [25l [3 El H] . Recently we have introduced a binary test for 
chaos, the 0-1 test, designed for the analysis of deterministic dynamical systems |8i|9j. 
The test distinguishes between regular and chaotic dynamics for a deterministic sys- 
tem. The nature of the dynamical system is irrelevant for the implementation of the 
test; it is applicable to data generated from maps, ordinary differential equations and 
partial differential equations. The test has been applied to noisy numerical data [9J, 
experimental data [5] , quasiperiodically forced systems and strange nonchaotic attrac- 
tors [3], Hamiltonian systems [26], nonsmooth systems [15] and fluid dynamics [T8] . 

The usual test of whether a deterministic dynamical system is chaotic or non- 
chaotic involves the calculation of the maximal Lyapunov exponent A [13] . A positive 
maximal Lyapunov exponent indicates chaos: if A > 0, then nearby trajectories 
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separate exponentially and if A < 0, then nearby trajectories remain in a close neigh- 
bourhood of each other. This approach has been widely used for dynamical systems 
whose equations are known. If the equations are not known or one wishes to examine 
experimental data, then A may be estimated using the phase space reconstruction 
method of Takens [27j, by approximating the linearisation of the evolution operator 
[21!, or by the "direct method" [23]. 

In contrast our test does not depend on phase space reconstruction but rather 
works directly with the time series given. The main advantages of our test are (i) 
it is binary (minimizing issues of distinguishing small positive numbers from zero), 
(ii) the nature of the vector field as well as its dimensionality does not pose practical 
hmitations, and (iii) it does not suffer from the difficulties associated with phase space 
reconstruction [T3] . 

In this paper, we describe in detail how to implement the 0-1 test for chaos. In 
addition, we carry out modifications to the test that greatly improve the previous 
versions in [H [9] . 

Throughout the paper, we use the logistic map to illustrate our claims, with 
the exception of Section [H] where we use the Lorenz attractor as an example of a 
continuous time system. The reader can verify that our results apply equally well to 
other systems, including those considered in our previous papers [H [9] . 

1.1 Recipe for the 0—1 test 

We briefiy review how the test is implemented. Given an observation for j = 
1, . . . , N we perform the following sequence of steps: 

1. For c G (0, vr), we compute the translation variables 



for n = 1,2, . . . , N . Typical plots of p and q for regular and chaotic dynamics 
are given in Fig. [T] 

2. The diffusive (or non-diffusive) behaviour of Pc and Qc can be investigated by 
analyzing the mean square displacement Mc{n). The theory behind our test 
assures that if the dynamics is regular then the mean square displacement is 
a bounded function in time, whereas if the dynamics is chaotic then the mean 
square displacement scales linearly with time. In Section [2] we look at expres- 
sions for the mean square displacement and describe how one may use analytical 
expressions derived in ^llj to conveniently modify the expression for the mean 
square displacement. 
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Figure 1: Plot of p versus q for the logistic map x„+i = /ix„(l — Left: Regular 
dynamics at fi = 3.55; Right: Chaotic dynamics at /x = 3.9. We used 5000 data 
points. 

3. We then compute the asymptotic growth rate Kc of the mean square displace- 
ment. Methods for the most effective estimation of this are discussed in Sec- 



4. Steps 1-3 are performed for Nc values of c chosen randomly in the interval 
(0, tt). In practice, Nc = 100 is sufficient. The choice of c is discussed further in 
Section m We then compute the median of these Nc values of Kc to compute the 
final result K = median(ii'c). Our test states that a value of K ~ indicates 
regular dynamics, and K ^ 1 indicates chaotic dynamics. 

In this paper, we explore practical issues arising in the implementation of the above 
algorithm. Various issues associated with steps 2-4 are discussed in Sections [2H1] 
respectively. In Section [5] we examine finite data size effects. In particular we look at 
weak chaos. In SectionEJwe consider continuous time systems where oversampled data 
can lead to small values of K despite an underlying chaotic dynamics. In Section [7] 
we investigate the issue of measurement noise. 

Remark 1.1 In the first version of our test, introduced in [8], we defined pdn) and 
qdn) by iterating the extended system 



tionO 



Pc{n + 1) 
qc{n + 1) 
Mn + 1) 



Pc{n) + (pin) cosi^cin)) 
qc{n) + (j){n) sm{dc{ri)) 
^c{n) + c + a (pin) . 
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The current version of the test corresponds to the case a = 0. As shown in [5], the 
test with a = is less sensitive to measurement noise. 

Remark 1.2 It can be rigorously shown that (i) pdn) and qdn) are bounded if the 
underlying dynamics is regular, i.e. periodic or quasiperiodic and (ii) Pc{n) and qdn) 
behave asymptotically like Brownian motion for large classes of chaotic dynamical 
systems. 

In [8] we used results of [6], |20l |2T] to prove this for the case a 7^ in Remark 11.11 
In the case a = these results are not applicable; nevertheless in [11] we cover the 
case a = under even weaker assumptions on the underlying dynamics. 

Remark 1.3 In [HI IE] it was shown that the results are valid for almost all observ- 
ables (p. Of course, the choice of the observable (p influences the rate of convergence 
(but not the limiting value K = or K = 1). From a practical point of view we 
found that changing the observable does not greatly alter the computed value of K. 



2 Computation of the mean square displacement 

For a given time series 0(j) with j = 1, . . . , A^, we compute the mean square displace- 
ment of the translation variables pd^) and qd^) defined in (II. ip for several values of 
c G (0,7r). The mean square displacement is defined as 

1 ^ 

Mdn) = lim - y^lpdj +n)- pdj)? + M + - qdj)? ■ (2.1) 

Note that this definition requires n <^ N. In [9] we calculated the mean square 
displacement using directly the definition (12.11) . The limit is assured by calculating 
Mdn) only for n < ricut where Ucut N. In practice we find that ricut = N/10 yields 
good results. 

The test for chaos is based on the growth rate of Mdn) as a function of n. In the 
following, we use analytical expressions derived in [11] to formulate a modified mean 
square displacement Ddn) which exhibits the same asymptotic growth as Mdn) but 
with better convergence properties. 

Under mild assumptions on the underlying dynamical system, described in Re- 
mark [2?T] below, for each c G (0,7r), 

Mdn) = V{c) n + Ksc(c, n) + e(c, n) , (2.2) 

where e(c, n)/n^Oasn^oo uniformly in c G (0, tt) and 

Ksc(c,n) = {E4>Y- . 

1 — cos c 
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The expectation E(j) is given by 

1 ^ 
i=i 

The form (12. 2p suggests an improvement for the test: We can subtract the exphcit 
term K)sc(c, n) from the mean square displacement and introduce 

D,H = M,(n)-Ksc(c,n). (2.3) 

Note that the asymptotic growth rates of Mc{n) and Dc{n) are the same. 



7.5- 




Figure 2: Plot of mean square displacement versus n for the logistic map with 
/i = 3.91 corresponding to chaotic dynamics. The oscillating (green) curve is the 
original mean square displacement Mc(n) as defined in (12. ip : the straighter (red) 
curve is the modified mean square displacement Dc{n) as defined in (12. 3p . We used 
2000 data points and computed Mc(n) and Dcin) for n = 1, . . . , 200 and c = 1.0. 

In Fig. [2] we show the two mean square displacements Mc{n) and Dc{n) for the 
logistic map Xn+i = /ia;„(l — x„) with /i = 3.91 (which corresponds to chaotic dynam- 
ics) and an arbitrary value of c = 1.0. Evidently, the subtraction of the oscillatory 
term K)sc(c, n) regularizes the linear behaviour of Mc{n). This allows a much better 
determination of the asymptotic growth rate Kc- 

Remark 2.1 The autocorrelation function for the observation 0(j) is given by 

pik) = Ei^imk + 1)) - {E^Y, k = 0,1,2... 

Provided the autocorrelations are absolutely summable (that is, YlT=o\p(^)\ ^ ^) 
then equation (12. 2p is valid, and moreover, the error term e(c, n) decays uniformly in 



5 



c G (0,7r) (see for example [H]). It is for this reason that the test based on Dc{n) 
greatly outperforms the test based on Mc{n). 

Furthermore, the absolute summability condition guarantees pjj that 

V{c) = V e^'^p{\k\) = hm -e| Ve^^-^0(j) (2.4) 

^ n-^oo n I ^ 



j=0 



for all c G (0, 2tt). This result follows from the Birkhoff ergodic theorem, the Wiener- 
Khintchine theorem, and standard calculations. In particular, the slope V{c) of the 
mean square displacement is identified with the power spectrum. 

For nonmixing systems, the error term e(c, n) no longer decays to zero and there 
are further oscillatory terms in addition to Vosc{c,n). Nevertheless, the identifica- 
tion fl2.4p remains valid for nonmixing systems under very weak conditions [T^ . 

More importantly from the point of view of the test for chaos, working with Dc{n) 
remains highly advantageous even for nonmixing systems. This is illustrated for the 
logistic map in Fig. H] later in this paper. 



3 Computation of Kc 

Having calculated the modified mean square displacement D c{n) foi n = 1,2,..., ncut, 
the next step is to estimate the asymptotic growth rate Kc. We have tried out 
two different methods: a regression method and a correlation method, described in 
subsections 13.11 and 13.21 below. 



3.1 Regression method 

The regression method consists of linear regression for the log-log plot of the mean 
square displacement. In [H] we used the original mean square displacement Mc{n), so 
the asymptotic growth rate Kc is given by the definition 

n^oo log n 

Numerically, Kc is determined by fitting a straight line to the graph of log Mc(n) 
versus logn through minimizing the absolute deviation [22] . 

In Section [2], we demonstrated the superiority of the modified mean square dis- 
placement Dc{n) when compared to Mc(n), so it is natural to apply the regression 
method to Ddn). Whereas Mc{n) is strictly positive, Ddn) may be negative due to 
the subtraction of the oscillatory term V^sc(c, n). Hence, we set 

Dc{n) = Dc{n) - min Dc{n) , 

n=l,...,ncut 
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and obtain the asymptotic growth rate 

n^oo log n 

Again, Kc can be determined numerically by regression (minimizing the absolute 
deviation) for the graph of log D(.{n) versus logn. 








-4 - 



2 4 6 

logn 

Figure 3: Plot of log Dc{n) as a function of logn for the logistic map at /i = 3.62. 
We used N = 2000 and calculated the mean square displacement up to ricut = A^/10. 

Remark 3.1 Minimizing the absolute deviation is preferable when compared to the 
usual least square method as the latter assigns a higher weight to outliers. Since the 
linear behaviour of the mean square displacement is only given asymptotically, one 
typically encounters outliers for small values of n. We find that it is usually sufficient 
to use the absolute deviation for estimating K^, and that it is not necessary to employ 
more complicated higher-order regression methods such as the method by Yohai [17]. 

The finite value of the o(n)-term e(c, n) in the definition of Mc{n) in (12.21) leads 
to a distortion for small values of n. In such situations, one typically observes a 
fiattening of the slope of log Mc(n) or log Dciji) as illustrated in Fig. [31 It is those 
values for small n of log Mc(n) (or log Ddn)) which would be overestimated in a least 
square fit. 

3.2 Correlation method 

We now present an alternative method for determining from the mean square 
displacement. (The method is described in terms of Dc{n), but we could use Mc(n) 
instead.) 



7 



Form the vectors ^=(1,2,..., ricut) and A = (Dc(l), Dc{2), . . . , Dc{ncut))- Given 
vectors x, y of length q, we define covariance and variance in the usual way: 



1 " 

cov{x,y) = - - x){y{j) - y), where 



1^ 



X = - 



i=i 



Q 

var(x) = cov(x, x) . 
Now define the correlation coefficient 

JC^co.xK.A)^ e 1-1.1] . 

A/var(4 jvar(A) 

This quantity measures the strength of the correlation of Dc{n) with linear growth. 
Again, it can be shown rigorously [llj that under weak conditions on the underlying 
dynamics (as described in Remark ll.2p we obtain Kc = for regular dynamics and 
Kc = 1 for chaotic djTiamics. 
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Figure 4: Plot of K versus /i for the logistic map with 3.5 < /i < 4 increased in 
increments of 0.001. We used 2000 data points. The darker (red) lines are obtained 
by using the original definition of the mean square displacement Mc{n) in fl2.ll) . The 
lighter (green) lines are obtained by using the modified mean square displacement 
Dc{n) in (12.31) . The resulting values of K are shown for the regression method (left) 
and the correlation method (right). The horizontal lines (blue and magenta) indicate 
the cases K = and K = 1. We used iY^ = 100 values of c. 

In practical terms, the correlation method greatly outperforms the regression 
method. This is evident from Fig. H] which compares the regression method and 
the correlation method (using both Mc{n) and Dc{n)) for the logistic map. 
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4 Choice of c and determination of K 



In Fig. [5] we show the asymptotic growth rate Kc as a function of c for regular and 
chaotic dynamics. In the case of periodic dynamics, most values of c yield i^"c = as 
expected, but there are isolated values of c for which Kc is large. (For the regression 
method, i^c ~ 2 at these resonant points.) These resonances are easily explained as 
follows: Equation (11. ip shows that if the Fourier decomposition of the observation 
contains a term proportional to exp(— ztu/c), then there is a resonance at c = u; where 
Pci^) ~ n, and hence Mc{n) ~ ra^, irrespective of whether the dynamics is regular 
or chaotic. For the plots in Fig. [5], we have calculated the asymptotic growth rate 
using both the regression method described in Section [XT] and the correlation method 
described in Section 13.21 
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Figure 5: Plot of K^ versus c for the logistic map calculated using the regression 
method (top) and correlation method (bottom). We used here N = 5000 data points, 
and 1000 equally spaced values for c. Left: /i = 3.55 corresponding to regular dynam- 
ics; Middle: fi = 3.9 corresponding to chaotic dynamics; Right: = 3.6 corresponding 
to chaotic but non-mixing dynamics. 



The occurrence of resonances for isolated values of c suggests using the median 
of the computed values of Kc- (We use the median rather than the mean, since the 
median is robust against outliers associated with resonances.) 

In Fig. [51:;, K^ is shown as a function of c for fi = 3.6 where the dynamics is chaotic 
but not mixing on the whole interval [0, 1]. The actual dynamics in the logistic map 
oscillates between two disjoint sets, each of which is mixing, and there is a resonance 



at c = TT. At resonance, pc 



n] 



n and Mc{n) ~ n as before. Close to resonance. 
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the p-q plot eventually behaves like Brownian motion, but in practice one sees only a 
small part of this motion and so Kf. ^ 0. 

Remark 4.1 Naturally, the choices of c are equally spaced in Fig. O whereas in 
applying the test (and throughout the paper with the exception of Fig. [5]) we choose 
randomly sampled values of c. 

To avoid that resonances distort the statistics, we further restrict the range of 
randomly sampled values for c to c G (7r/5,47r/5) for all our computations. The 
resonance at c = is inherent to our test, but it may leak through adjacent values of 
c as seen in Fig. [St. The further restriction to exclude vr is not necessary, but we found 
it helpful. (A typical route to chaos is the Feigenbaum route via period doubling. 
Here, the parameter ranges for fixed points and period two points are largest.) 

In Fig. [HI we show how the result for K depends on the number A^^^ of different 
values of c. Here we use the correlation version of the test to calculate Kc as described 
in Section 13. 2[ There is no measurable gain in increasing Nc from 100 to 1000 and 
we find that generally Nc = 100 different values of c is sufficient. 

5 Finite size problems 

There are three types of finite size effects. First, the time series needs to be long 
enough to explore and sample the relevant phase space area (i.e. the attractor). This 
is an inherent problem affecting all tests for chaos. Second, the definition of the mean 
square displacement involves a limit which requires n <^ N. Accordingly, we have 
chosen n < ricut = N/10. 

Third, the theory developed in [8], |TT] makes statements about the asymptotic be- 
haviour of Dc{n) (or Mc{n)) and as such requires ricut, and hence A^, to be sufficiently 
large. Finite size effect in this context means that for small n the asymptotic linear 
growth is not yet dominating, see Fig. [31 This finite size effect is explored in the 
remainder of this section. From now on, we work exclusively with the modified mean 
square displacement Dc{n) and the correlation method. 

In Fig. [7] we show how the value of K depends on the amount of data used. We 
can see clearly the convergence towards the asymptotic values K = and K = 1 for 
regular and chaotic underlying dynamics, respectively. (For values of /i corresponding 
to stronger chaotic dynamics well within the chaotic range, the convergence towards 
A' = 1 is even more rapid.) 

In the case of "weak chaos", close to the so called "edge of chaos", longer data 
sets are required to obtain K = 1. Weak chaos is characterized by a slow decay 
of correlations. This has consequences for the modified mean square displacement 
Dc{n) = V{c)n+o{n). For systems whose auto-correlation function is slowly decaying. 
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Figure 6: Plot of K versus for the logistic map using the correlation method, with 
3.5 < < 4 increased in increments of 0.001. We used 2000 data points. Upper left: 
= 1; Upper right: = 10; Lower left: = 100; Lower right: A^^ = 1000. 



it may be the case that the o{n) term dominates for the available data. We illustrate 
this problem in the context of the logistic map. The bifurcation parameter fi takes 
the value /x = /Xoo = 3.569945672 ... at the edge of chaos and for /i = /ioo + 0.001 one 
observes weak chaos. 

It has been erroneously claimed that our test cannot detect weak chaos, see [T2|[T0]. 
In fact, there are two methods whereby we can distinguish between regular dynamics 
and weak chaos: 

(i) By visual inspection of the plot in the p-q plane generated as in Fig. [HI (Note 
that for longer data sets the dynamics in the p-q plane in Fig. [Hb would look 
just like Fig. [Tb.) 

(ii) By looking at the dependence of as a function of A^. As illustrated in Fig. [9l 
we can distinguish weakly chaotic from regular dynamics even when the value 
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Figure 7: Plot of K versus the available amount of data for the logistic map. Left: 
II = 3.55 corresponding to regular dynamics; Right: /i = 3.6 corresponding to chaotic 
dynamics. 

of K is very small - note that K = 0.027 for = 2000 in the weakly chaotic 
case. 
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Figure 8: Plot of p versus q for the logistic map. Left: /x = Hoo] Right: n = /ioo+0.001. 
We used 5000 data points. 
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Figure 9: Plot of as a function of for the logistic map at /i = /ioo (left) and 
At = /Uoo + 0.001 (right). Although the value of K is small in both cases, the behaviour 
ofK function of N distinguishes the two cases. 
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6 Continuous time systems 



In the previous sections, the 0-1 test was formulated for discrete time systems. For 
continuous time series there is a well-known oversampling issue that must be 
addressed. In this section, we discuss this difficulty and how to overcome it. 

Given < ti < t2 < ^3 < ■ ■ ■ we obtain a discrete time series 0(^2), ^(^3), • • • 

to which the test for chaos may be applied as in previous sections. (The sequence 
tj, j > 1, should be chosen in a deterministic manner so that the time series (f){tj) 
is deterministic.) One method of choosing the tj is as the intersection times with a 
cross-section, so the time series (f){tj) corresponds to observing a Poincare map. In 
this situation, there are no issues with oversampling. 

A second, perhaps more usual, approach is to take tj = jTg where > is the 
sampling time. The time series = (pijTg) corresponds to observing the "time-Ts" 
map associated with the underlying continuous time system. If r<j is too small, then 
the system is oversampled and this often leads to incorrect results. To illustrate the 
issue of oversampling we study the 3-dimensional Lorenz system 

X = 10{y — x) 

y = 30 x — y — xz (6.1) 
z = xy-^z , 

which exhibits robust chaos. We have integrated this system with a time step of 
At = 0.001 and recorded 100,000 data points (ie. 100 time units). 

Fig. [To] shows an oversampled and a sufficiently coarsely sampled observable for 
the Lorenz system (16. ip . The finely sampled time series (r^ = 0.005) yields K ^ 
even for = 100, 000 whereas the coarsely sampled data (ts = 0.05) yields K ^ 1 
already for = 5, 000. 

A good choice of the sampling time can often be obtained by visual inspection 
as in Fig. [TOl A more refined method is to use the first minimum of the mutual 
information [H |13]. For the data depicted in Fig. [10] this method yields = 0.17 
(roughly a quarter of the oscillation period). Note however that in this particular 
instance the smaller sampling time = 0.05 already gives K ^ 1 and extracts a 
longer time series from the data in Fig. [10] In general, the optimal sampling time 
will depend on the dynamical system and the time series under consideration. We 
refer the reader to [13] for a discussion on optimal time delays in the context of phase 
space reconstruction. 

Although oversampling is a practical problem for data series of finite size, it should 
be emphasized that theoretically the test works for all sampling times Tg in the limit 
N ^ 00. 
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Figure 10: Plot of the observable (j){t) = x{t) for the Lorenz system (16. ip . The finely 
sampled data (red) are sampled at = 0.005 time units. The coarsely sampled data 
(green filled circles) are sampled at = 0.05 time units. 



6.1 Oversampling and power spectra 

For continuous time systems, the mean square displacement is defined as 

M,(t) = hm /" M + r) - p,{r) f + {q,{t + r) - q,{T)f dr . 
For a time series sampled with sample time this can be approximated by 



1 ^ 

Me(n) = lim - ^ ([per.(j + u) - Pcr^j)]' + [QctAJ + ^) " Qcr^j)?) 
iV — >oo I\ ' 



Similarly the power spectrum for the time-continuous case discretizes to 

n-l 



S{u) = lim iE|^e^-^^V(j) Vf, (6.2 



j=Q 



where z/^ = l/r^ is the sample frequency. The power spectrum consists of discrete 
peaks if the underlying system is regular, and is nowhere zero for a large class of 
chaotic systems fi9\. However, for chaotic systems the power spectrum decays for 
large frequencies u, and so for frequencies larger than some z/max the power spectrum 
is zero for all practical purposes. 
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Comparing fl6.2p with the power spectrum fl2.4l) for discrete-time data, we identify 



C = 27r — , V e [0, Z/max] • 

Samphng at the Nyquist rate with = 2z/max yields c G (0, n) as before. However, 
oversamphng at a higher frequency Us > u*, restricts the effective choices of c to 
c G (0, c*) where c* = ^vr < tt. There is now a positive probabihty that the test 
for chaos will incorrectly yield K = since it is possible that more than half of the 
randomly chosen values of c G (0, vr) will lie in (c*, vr). 

We illustrate the previous argument using the Lorenz system (16.11) sampled with 
Ts ranging from = At up to = 300At. In Fig. [11] the median of the asymptotic 
growth rate K is shown as a function of the sample time. For data that is too 
finely sampled, we obtain K = although the dynamics is actually chaotic. Fig. [12] 
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Figure 11: Plot of i^' as a function of the sample time Tg for the Lorenz system (16. ip . 
The sample time is measured in units of At = 0.001. 

illustrates how the range of effective values of c depends on the sampling time r^. 

7 Noise contaminated data 

Real- world data is invariably contaminated with noise. Any method for distinguishing 
regular from chaotic dynamics can only succeed if the noise-level is sufficiently small. 
There are various standard noise reduction techniques [13j that may be applied in 
advance of applying any given test for chaos. In addition, the test itself may be 
modified. Below we indicate a modification of the 0-1 test for chaos that makes it 
more robust to the presence of noise. 

In [S] we introduced a version of the 0-1 test that works well for data contaminated 
with measurement noise. (This is the test as presented in Sections [2] and [31 but using 
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Figure 12: Plot of as a function of the frequency c for the Lorenz system (16. ip . 
From left to right we used r, = 5At, r, = lOAt, = 20At, r, = 30At, t, = 50At, 
Ts = 70At. The linear scaling of the range of c for which Kc ^ 1 is evident in the 
relative spacing of the respective lines. 

Mc{n) instead of Dc{n) and using the regression method instead of the correlation 
method.) We showed that the 0-1 test clearly outperforms tangent space methods 
and compares favourably to "direct methods" based on phase space reconstruction. 
The improvements in this paper have made our test extremely sensitive to weak 
chaos. However, an unavoidable consequence is an increased sensitivity also to noise 
(see Fig. [13] below). 

It turns out that the success of the version of the test in [9] is due to the oscillatory 
term V^sc(c, n) = (Ecj))'^ \"lTosT ^^^^ subtracted in Section [2] to define the modified 
mean-square-displacement Dc{n) = Mc{n) — Vosc{c,n). This term desensitizes the 
test and damps the ability to detect slow growth of the mean-square-displacement for 
time-series data of moderate length. Instead of reintroducing this term we adopt a 
more flexible approach, defining 



(The frequency \/2 was chosen arbitrarily.) For a large, we expect K = 0. The 
amplitude a of the term Vdamp(^) controls the sensitivity of the test to weak noise 
and simultaneously to weak chaos. This trade-off is unavoidable in any test for chaos. 

As an illustration, we consider the logistic map with measurement noise. Take as 
observable (f){n) = x„ and write 



where rjn are i.i.d. random variables drawn from a uniform distribution on [—1, 1] and 
e is the noise-level in percent. Fig. US] shows how the undamped version of the test 



D*^{n) = Dc{n) + aVdampW, KiampH = {E(t)f sin{V2n) . 



^{n) = 0(n)(l + —Vn) 
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in this paper copes with a noise level of 10% and the improvement that is obtained 
by using the damped mean-square-displacement D*{n). We obtain similar results for 
normally distributed noise. 
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0.4 
0.2 ^ 





LJ 
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Figure 13: Plot of K versus /i for the logistic map increased in increments of 0.001. 
The darker (red) lines were computed using clean data. The lighter (green) lines 
were computed after addition of 10% uniformly distributed measurement noise. Both 
lines were computed using the undamped mean-square-displacement Ddn). Left: 
= 1000 using the undamped mean-square-displacement Dc{n), Right: = 5000 
using Dc{n). Bottom: N = 5000 using the damped mean-square-displacement D*{n) 
with a = 2.5. 



Remark 7.1 Under the assumption that the noise is diffusive and not correlated 
with the dynamics, the mean square displacement for data contaminated with mea- 
surement noise may be written as 

Dc{n) = (Vdyn(c) + Koisc(c))n + o{n) , 
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where for a given value of c, Vdyn(c) is the variance associated with the deterministic 
dynamics and Vnoise(c) the variance associated with the measurement noise. Consider 
an ideahzed situation where the value of V^oise(c) is roughly constant as a parameter 
A is varied. Suppose further that the dynamics is known to be regular at A = Aq. 
Then we may estimate Vnoise(c) by making a gauge- measurement at A = Aq, applying 
the correlation method to Dc{n) — Vn. The unique value V = K(Ao) which yields 
Kc = is our estimate for l^oise(c). For other values of A we may now apply the 
correlation method to Dc{n) — Vc{Xo)n. 



8 Discussion 



We have presented a guide for the implementation of the 0-1 test for chaos. At the 
same time, we have introduced an improved version of the test which uses analytical 
expressions derived in [11]. Issues such as oversampling for continuous-time data and 
the presence of noise have been discussed. We hope that this guide will be helpful for 
scientists who would like to use the test. 

There are numerous methods in the literature for distinguishing between deter- 
ministic and chaotic dynamics. In our previous papers [8l [9], we made a careful 
comparison of the 0-1 test with methods for computing the maximal Lyapunov ex- 
ponent. Another method is to use the power spectrum for which there are efficient 
computational techniques. It should be pointed out however that these techniques 
generally rely on the Wiener-Khintchine theorem which assumes summable decay of 
correlations and hence excludes periodic and quasiperiodic dynamics. Hence to use 
power spectra as a test for chaos, it seems necessary to avoid the Wiener-Khintchine 

n-l ■• ^ 

theorem and to work directly with the expression lim„^oo X]j=o ^^^'^'PH) ■ Equa- 
tion ( 12. 4p shows the relationship between the 0-1 test and power spectra, and our 
test can be viewed as a way of condensing the information relevant for chaoticity or 
regularity contained in the power spectrum into a single binary number. 
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